Optical control and entanglement of atomic Schrodinger fields 
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We develop a fully quantized model of a Bose-Einstein condensate driven by a far off-resonant 
pump laser which interacts with a single mode of an optical ring cavity. In the linear regime, the cavity 
mode exhibits spontaneous exponential gain correlated with the appearance of two atomic field side- 
modes. These side-modes and the cavity field are generated in a highly entangled state, characterized 
by thermal intensity fluctuations in the individual modes, but with two-mode correlation functions 
which violate certain classical inequalities. By injecting an initial coherent field into the optical cavity 
one can significantly decrease the intensity fluctuations at the expense of reducing the correlations, 
thus allowing for optical control over the quantum statistical properties of matter waves. 

PACS numbers: 03.75-b,03.75.fi,42.50.Vk,42.55-f 
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The recent demonstration of Bose-Einstein condensa- 
tion in low-density alkali vapors [QJ^] opens up a new 
paradigm in AMO physics. It is now possible to gen- 
erate macroscopic atomic fields whose quantum statisti- 
cal properties can in principle be manipulated and con- 
trolled, very much like those of quantum optical fields. 
One important consideration is to determine to which 
extent the quantum state of a many-particle atomic field 
can be optically manipulated. In the single-particle case, 
the answer to this problem is known to a large extent. 
This is the domain of atom optics where a number 
of optical elements for matter waves have now been de- 
veloped, including gratings, mirrors, interferometers, res- 
onators, etc. But these optical elements manipulate just 
the atomic field "density" , or at most first-order coher- 
ence properties. However, Schrodinger fields possess a 
wealth of further properties past their first-order coher- 
ence, including atom statistics, density correlation func- 
tions, etc. In analogy to the optical case, one can there- 
fore think of "quantum atom optics" as that extension of 
atom optics where the quantum state of a many-particle 
matter-wave field is being controlled, characterized and 
used in novel applications. 

This Letter presents an analysis of a system where an 
optical field is used to manipulate the quantum state 
of a matter-wave field. Using a geometry identical to 
that used in the collective atomic recoil laser and 
very similar to those used for atom interferometry 0, 
and recoil-induced resonances |]J9], we consider a Bose- 
Einstein condensate driven by a far off-resonant pump 
laser and coupled to a single mode of an optical ring- 
cavity. This results in gain in the cavity mode, as well 
as the generation of momentum side- modes of the BEC, 
which are assumed to be unbound by any magnetic or 
optical trap. The quantum statistical properties of the 
side-modes can be strongly manipulated by varying the 
initial state of the optical cavity mode. In addition, a 
strong quantum mechanical entanglement can develop 
between the optical and matter-wave fields, as well as 



between matter-wave side-modes. The experimental re- 
alization of this system is currently feasible, in view of 
recent experiments on the diffraction of condensates by 
the NIST group (hJ. 

We consider an ultracold sample of bosonic atoms 
driven by a strong classical "pump" and a counterpropa- 
gating weak quantized "probe" optical field, both being 
far off-resonant from any electronic transition. Under 
these conditions the internal atomic degrees of freedom 
can be adiabatically eliminated and the matter- wave field 
is effectively scalar. The combined Hamiltonian for the 
atomic and probe fields is 
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Here, fio is the Rabi frequency of the pump laser of fre- 
quency wo an d momentum ko, A is the annihilation op- 
erator of the probe field of frequency u> and momentum 
k, satisfying [A, A*] = 1, and c(q) is the annihilation op- 
erator for a ground state atom of momentum q, satisfy- 
ing [c(q), c T (q')] = £q, q '. In addition, A is the detuning 
between the pump frequency and the upper electronic 
level closest to resonance, and g = d[ck/(2eoLS)] 1 ^ 2 
is the atom-probe coupling constant. Here d is the 
atomic dipole moment, L the length of the ring resonator 
sustaining the probe mode, and S the cross-section of 
that mode in the region of the atomic sample. Finally, 
K = ko — k is the atomic recoil momentum resulting 
from the absorption of a pump photon followed by the 
emission of a probe photon. 

The first two terms in Eq. (Q) are the free Hamilto- 
nians of the atomic and probe fields, respectively. The 
remaining terms correspond to the various processes by 
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which an atom undergoes a virtual transition under the 
influence of the optical fields. The first such term in- 
volves the exchange of a photon between the pump and 
probe fields, e.g. stimulated absorption of a pump pho- 
ton followed by stimulated emission of a probe photon, or 
vice versa. In coordinate representation, this term would 
take the form of the familiar periodic optical potential 
generated by the counterpropagating pump and probe 
lasers fields. The last two terms in Eq. (|l]) correspond 
to processes where a photon is first absorbed and then 
reemitted into the same field. These transitions are re- 
coilless, but contribute a cross-phase modulation between 
the atomic and optical fields. 

Assuming that the initial momentum width of the con- 
densate is small compared to the recoil momentum K, it 
is reasonable to treat it as a single mode atomic field of 
momentum q = 0. We furthermore restrict our discus- 
sion to the case TCT C , where T c is the critical temper- 
ature, and assume a large condensate for which the bare 
mode q = can then be described to a good approxima- 
tion as a c-number, c(0) — > V^/V exp(i|f2 | 2 ^/4A), where 
N is the mean number of atoms in the condensate. This 
approximation neglects both the depletion which occurs 
as atoms are transferred into the side-modes q ^ and 
the cross-phase modulation between the condensate and 
the probe field, thus it is valid for times short enough 
that E q#0 (ct(q)c(q)) « N, and (iU) « |fi | 2 /4| ff | 2 . 
This is the matter-wave optics analog of the familiar 
classical and undepleted pump approximation of nonlin- 
ear optics. Hence we describe the optical and matter- 
wave fields on equal footings, treating all strongly popu- 
lated modes classically and all weakly populated modes 
quantum- mechanically. 

Once we have replaced the condensate mode with its 
c-number counterpart, we then neglect all terms in the 
Hamiltonian ([!]) involving the product of three or more 
weakly populated field modes. This is a direct conse- 
quence of Bose enhancement, which strongly strengthens 
the interactions involving the central mode q = mode 
relative to those involving only the side- modes, and leads 
us to the effective Hamiltonian 
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where uj r — hK 2 /2m, and we have introduced the 
slowly varying operators c± = exp(i|£7 | 2 £/4A)c(±K) 
and fi = — i(g*^A/|g||Slo||A|) exp(itOot)A. The system 
is fully characterized by the effective coupling constant 
X = | .9 1 |^o I V^V /8uvA and the dimensionless pump-probe 
detuning 8 = (u>o — uj)/uj r . 

The Hamiltonian (Q) describes three coupled field 
modes: the optical probe and two atomic condensate 
side-modes with wavenumbers ±K. The presence of 
terms such as a)c*_ in Eq. (j^) is a direct consequence 





of momentum conservation, and does not result from 
keeping antirotating terms in the Hamiltonian. It imme- 
diately brings to mind the optical parametric amplifier 
jll| , a device known to generate highly non-classical op- 
tical fields exhibiting two-mode intensity correlations and 
squeezing, and which has been extensively employed in 
the creation of entangled photon pairs for fundamental 
studies of quantum mechanics, quantum cryptography 
and quantum computing. A novel aspect of the present 
system is that it offers a way to achieve quantum en- 
tanglement not just between optical fields, but rather 
between macroscopically populated atomic Schrodinger 
and optical fields. 

The dynamics of the system can be determined by solv- 
ing the three coupled-mode equations 
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where t = uj r t, and we have introduced d a (r) = a(r), 
d-(r) = c_(r), and d+(r) = c + (r) for future notational 
compactness. 

An analytic solution can be constructed explicitly from 
the eigenvalues {Xj} and eigenvectors {vj} of the matrix 
on the right-hand side of (||). The eigenvalues {Xj} have 
been studied in detail in Ref. in the context of the 
theory of the Collective Atom Recoil Laser (CARL). It 
was shown that provided the system parameters x an d $ 
satisfy certain threshold conditions, they take the form 
Ai = wi, X2 = A3 = fi+ir, where u>i, f2, and T are all real 
quantities. Hence we see that after an initial transient, 
the solution grows exponentially in time at the rate T. 
This regime of exponential growth is familiar from the 
physics of the free-electron laser and of the CARL, where 
it is usually studied at high temperatures. The explicit 
form of the eigenvalues and eigenvectors is not required 
for the current analysis and will be presented elsewhere. 
For our present purposes, it is sufficient to know that for 
a given set of parameters Xi $ they are simply constants. 

The solution of Eqs. (||) is 



where the coefficients (r) are given by 
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Here V -1 is the inverse of the matrix V whose elements 
Vij are given by the ith component of the eigenvector v,- , 
and Qj = v^v^ 1 . The approximate equality in Eq. (||) is 
valid for times long enough that one can neglect all but 
the exponentially growing terms, henceforth referred to 
as the exponential growth regime. 
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This exponential growth of the system can be triggered 
either from vacuum fluctuations, as we discuss in more 
detail shortly, or by a weak injected probe signal. We 
investigate both situations by assuming that the probe 
field is initially in the coherent state a, the vacuum state 
corresponding to a = 0. The condensate side-modes, in 
contrast, are always taken to be in the vacuum state at 
t = 0, so that the initial state of the system is \a, 0, 0). 

The expectation values (di) of the three coupled modes 
are readily found to be 



(di(r)} = au ia (T) w a( ia e 



(T+i(l)i 
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where the approximate result is for the exponential 
growth regime. As expected, in the absence of an injected 
signal the mean fields remain zero, but the injected probe 
breaks the symmetry of the system and lead to nonzero 
expectation values. 

Decomposing the expectation values (di) in terms of an 
amplitude and phase as (di) = £i(r) exp[i^»j(r)], we find 
that in the exponential growth regime their uncertainties 
obey 
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where the fluctuation function f(x, 5) — \ v 3-/ v 3 a \ nas a 
relatively simple dependence on the control parameters 
X and 6. Specifically, for a given f(x, 5) is approxi- 
mately unity at the <5 which maximizes the growth rate 
r, and increases steadily away from this value. Clearly, 
Eq. (R) holds only in the case of an injected probe signal, 
a 7^ 0. In that case, the phase uncertainties of all three 
mean fields approach the same limiting value for large r, 
and this value approaches zero as a becomes very large, 
i.e. for large enough a all three modes are effectively in 
coherent states. We note that for large a the system is 
essentially equivalent to Kapitza-Dirac atomic diffraction 
of a condensate by a standing wave flOl , which we now 
see produces atomic side-modes in coherent states. 

We now investigate the mean intensities = 
(d|(r)di(r)) — The 5-function accounts for the fact 
that d%(T) in Eq. (|J) is a creation rather than an annihila- 
tion operator, thus guaranteeing that the initial intensity 
of the"-" side-mode vanishes. These intensities are given 
explicitly by 



h(T) = \a\ 2 \u ia (T)\' + k_(7 
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In the exponential growth regime they reduce to 

Ii(r)^(\a\ 2 \C ia \ 2 + \Q-\ 2 )e 2 ^. (9) 

They have a stimulated component, proportional to |a| 2 , 
and a spontaneous component, which is present even 
when all three field modes begin in the vacuum state. 
The stimulated component is simply the squared ampli- 
tude of the mean field, while the spontaneous component 



has no mean field, as it originates from the amplification 
of vacuum fluctuations in the atomic bunching. 

To help understand this in more detail, we in- 
troduce the atomic "bunching operator" B = 
(1/N) J2j exp(iKzj), where ij is the position operator 
of the jth atom. If the atoms in the sample are evenly 
distributed in space then (B) =0. At the opposite ex- 
treme, if all the atoms are localized on a array of pe- 
riod 2ir/K, then \(B)\ = 1. Second-quantizing B and 
linearizing the result by treating the q = mode as a 
c-number and keeping as in the derivation of Eq. (j^) 
only the lowest order terms in the side-mode operators, 
we can reexpress B in terms of the atomic field operators 
as B = (l/y/~N)(cL + c+). It is immediately apparent 
from that definition that (5(0)) = for our intial state 
| a, 0, 0). However, the fluctuations (B 2 (0)) are nonzero, 
due to the fact that (c_(0)cl(0)) = 1. It is precisely 
this expectation value which leads to the spontaneous 
intensity component, which can therefore be attributed 
to vacuum fluctuations in the initial atomic bunching. 
These fluctuations play a role similar to that of vacuum 
fluctuations in spontaneous emission. 

In addition to the side-mode intensity, it is instruc- 
tive to also study their equal-time intensity correlation 
functions, as they provide useful information about their 
quantum or classical nature. For the probe mode, we 
have 



fflt(T)fflt(T)a(r)a(r)) 



(10) 



(2) (2) 

The side-mode correlation functions g_ (r) and g\ (r) 
are defined likewise but with a(r) replaced by c_(r) and 
c+(t) respectively. These correlation functions are given 
explicitly as 



a 2 \u la (T)\ 4 
If(r) 



N 2 + /(x,<5)F 



(ii) 



As before, the approximate result applies to the expo- 
nential growth regime, where the intensity correlation 
functions become constant in time and the same for each 
mode. In the case where the system builds up from noise 
(\a\ 2 = 0), we have then = 2, the signature of a 
thermal or chaotic field. As the injected signal strength 
is increased, however, g^ — > 1, which is characteristic 
of a Glauber coherent field with Poissonian excitation 
statistics. Note the important point that the state of the 
side-modes can be continuously varied from thermal to 
coherent by varying the strength of the injected probe sig- 
nal and/or the system parameters x an d S. This demon- 
strates that the state of coherence of the matter-wave 
field can be directly controled by an optical field, a new 
form of coherent control. 

We have mentioned the analogy between the problem 
at hand and the parametric oscillator. It is the tool of 
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choice for generating entangled quantum optical states. 
We now investigate if similar entanglements can be ob- 
tained here. We proceed by investigating the equal-time 
two-mode intensity cross-correlations, which are a mea- 
sure of the degree of entanglement between the modes of 
the system. For example, the intensity cross-correlation 



function g a _(r) is defined as 

(at(r)a(r)cl(r)c_(T)) 



(2) 
7a- 



(at(r)a(r))(cL(r)c_(r)) 



(12) 



Other intensity cross-correlation functions such as g a \ (r) 

( 2) 

and g_+(r) are defined similarly. 

For classical fields, there is an upper limit to the 
second-order equal-time correlation function. It is given 
by the Cauchy-Schwartz inequality |I1J 
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(r)< gf'{r) g)"(r) 
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Quantum mechanical fields, however, can violate this in- 
equality and are instead constrained by the inequality 

m 



9$\t) < 
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Ii(r) 



1/2 



(14) 



which reduces to the classical result in the limit of large 
intensities. 

We focus our attention on the spontaneous case a = 
0, where the single-mode intensity correlation functions 
are given by g\ (r) = 2. In this case, the equal-time 
intensity cross-correlation functions are found to be 
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(15) 
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From Eq. (15) we see that both g a _(r) and g_+(r) vio- 

(2) 

late the Cauchy-Schwartz inequality, while g a _l_(r) is con- 
sistent with classical cross-correlations. Furthermore, the 
explicit evaluation of the £ij 's shows that I + (r) -C I a (t) , 
which implies that g a _ (r) is very close to the maximum 
violation of the classical inequality consistent with quan- 

(2) 

turn mechanics, whereas for g_i_(T) the violation is not 
as strong. In the two-mode parametric amplifier, the 
two-mode correlation function shows the maximum vio- 
lation of the Cauchy-Schwartz inequality consistent with 
quantum mechanics. In the three-mode system, however, 
the two-mode cross-correlation functions involve a trace 
over the third mode, hence it is not surprising that the 
two-mode correlations are not maximized. 

If we now allow for an injected coherent probe field 
(a =/= 0), we must first note that the intensities are in- 
creased by approximately |a| 2 , which means that the 



time scale on which the classical and quantum upper lim- 
its ( fl3| ) and ( |i~4| ) converge is reduced by l/|a| 2 , making an 
experimental confirmation of quantum correlations more 
difficult. In addition, whereas for the spontaneous case 

(2) 

a = 0, numerics show the cross-correlation g y a _ follows 
almost exactly the quantum upper limit ( |l4]) for all f > 0, 
for a 7^ 0, it lies somewhere in between the quantum jl^ ) 
and classical (13) limits. As a is increased, it falls ever 
closer to the classical upper limit, so that in the limit 
of very large a, the fields exhibit only classical cross- 
correlations. 

In summary, we have discussed how the quantum state 
of momentum side-modes of a condensate can be varied 
continuously between two distinct limits by specifying 
the initial state of an optical cavity mode. When it begins 
in the vacuum state, the side-mode and the cavity mode 
fields develop with zero mean fields, thermal intensity 
fluctuations, and strong quantum correlations between 
the modes. In contrast, when it is prepared in strong 
coherent state, we approach a "classical" limit in which 
the fields develop with non-zero mean fields having well 
defined phases, intensity fluctuations indicating a coher- 
ent state, and exhibiting classical correlations only. Con- 
densate side-modes have recently been realized in an ex- 
periment at NIST, where a condensate was subjected to 
Kapitza-Dirac diffraction by a standing-wave laser field. 
fiof In order to observe the effects predicted here, this 
field would need to be replaced by a combination of a 
strong pump laser and a weak counter-propagating probe 
sustained by an optical cavity, which does not appear to 
present any major difficulty. 
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